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Abstract: We investigate whether coherent flavor conversion of neutrinos in a neu- 
trino background is substantially modified by many-body effects, with respect to the 
conventional one-particle effective description. We study the evolution of a system of 
interacting neutrino plane waves in a box. Using its equivalence to a system of spins, 
we determine the character of its behavior completely analytically. We find that, if the 
neutrinos are initially in flavor eigenstates, no coherent flavor conversion is realized, in 
agreement with the effective one-particle description. This result does not depend on 
the size of the neutrino wavepackets and therefore has a general character. The validity 
of the several important applications of the one-particle formalism is thus confirmed. 
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1. Introduction 

The flavor composition of a neutrino system may be modified as a result of the interac- 
tions between the neutrinos and a background medium. The efficiency of this process 
depends, in addition to the properties of the background itself (composition, density, 
temperature, etc.), on the coherent (or incoherent) character of the interaction. 

Here we study the case of coherent neutrino conversion in a background of neutri- 
nos. The interaction of a neutrino beam with a neutrino background has the peculiar 
feature that, due to momentum exchange, a neutrino from the background may be 
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scattered into the beam and vice versa. This differs from the case of scattering on elec- 
trons and/or nucleons, and implies that the effect of coherent scattering on neutrinos 
can not be treated in analogy with that on ordinary matter. 

In the literature, the effects of the neutrino-neutrino coherent scattering on neutrino 
conversion have traditionally been treated following the pioneering work by Pantaleone 
[|IJ, EJ, Sigl & Raffelt and McKellar & Thomson M. In these papers, a single- 
particle evolution equation for each neutrino is derived, which has the form of the usual 
oscillation equation for neutrinos in matter with an extra term designed to capture 
the effect of neutrino-neutrino scattering. According to this equation, the effect of 
the interaction with the neutrino background is significant if the neutrino-neutrino 
potential is at least comparable to the potential due to matter (zero-temperature and 
thermal terms) and/or to the vacuum oscillation terms: 

{\n e - - n e +\ 
(n e - + n e +)El/m 2 w (1.1) 
Am 2 /(G F E U ) 

Here n v (rip) is the (anti)neutrino number density, n e - (n e +) is the electron (positron) 
number density, E v is the neutrino energy, m w is the W^-boson mass and Gf is the 
Fermi constant. 

The condition (|1 . 1|) may be satisfied in the early Universe close to the Big Bang 
Nucleosynthesis epoch (at temperature T ~ 0(1) MeV) or just outside a supernova 
core (at r ~ 0(100 km)). Indeed, in these environments the number density of neu- 
trinos is comparable to, or larger than, that of electrons and nucleons. At the same 
time, incoherent scattering is negligible and neutrinos stream freely in the medium, 
affected by coherent scattering (refraction) only. The evolution of the neutrino flavor 
composition in both cases has been extensively studied using the one-particle equa- 
tion. An important result [§, [| |7], ||, ^ is the strong bound on the chemical potential 
of non-electron relic neutrinos that follows from studying neutrino oscillations in the 
early Universe, with important implications on BBN. Other important applications in- 
clude hypothetical oscillations between active and sterile neutrino states that may have 
generated a lepton asymmetry in the early Universe [ID[ [II, 12, |T^, [14| and possible 



effects of a sterile neutrino on the r-processes in a supernova JER ITB, O, [L8|, [19] 



There is an important theoretical issue at the foundation of all of the mentioned 
analyses that, if unresolved, casts doubt on their validity. The issue is the existence and 
the range of validity of the description of neutrino self-refraction in terms of a set of 
one-particle equations. Generally speaking, the interactions between neutrinos in the 
ensemble may be expected to create quantum correlations, or, in other words "entan- 
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gled" states of many neutrinos which are not products of individual wavefunctions. A 
priori, a system in which such entanglement exists requires a many-particle description. 
The importance of this point was recognized early on [fj, but was not pursued further 
in the following years. 

Recently, the problem of possible existence and effects of quantum correlations 
has been investigated using two different approaches. The first, put forth in our recent 



paper ||20|| , is to understand the flavor evolution of a many-neutrinos system as an effect 
of the interference of many elementary neutrino-neutrino scattering events. Using this 
construction, we argued that, in the limit of infinite number of neutrinos, the many- 
body description factorizes into one-particle equations. These equations coincide with 
those given in the previous literature up to terms which change the overall phases 
of neutrino states. While such terms can be potentially important in non-oscillation 
phenomena, they do not affect the flavor evolution of the system. 

As a particular example, we considered an ensemble of neutrinos which are initially 
in flavor eigenstates. Following our analysis, we found that in this system, somewhat 
counterintuitively, no coherent flavor conversion takes place. This finding is in agree- 
ment with the standard one-particle description of |], H, H, EH- 

The second approach, introduced by Bell et al. , is to consider the evolution of 
a system of interacting neutrino plane waves in a box. The idea is to study the proper- 
ties of the many-body Hamiltonian of this system with the goal of determining the rate 
with which statistical equilibration is achieved. The analysis is applied to a system of 
neutrinos which are initially in flavor eigenstates. It is concluded that the time scale t eq 
of the flavor evolution of a given neutrino in this system is inversely proportional to the 
neutrino density and to the Fermi constant: t eq oc (nGV) -1 , which is characteristic of 
coherent flavor conversion. This result contrasts with the prediction of the one-particle 
description, and therefore has been considered as an indication of its breakdown and 
of the existence of new coherent effects due to entanglement. Clearly, this would have 
profound implications on the applications of neutrino-induced neutrino conversion to 
supernova physics and cosmology. 

Given the importance of the subject, in this paper we present a further study of 
neutrino flavor conversion in a neutrino background. The aim is to clarify the origin 
of the contrast between Ref. [[H]] and the other literature and give a definite answer to 
the question of the importance of entanglement and its effects. 

Since the suggested new effect comes from simultaneous interaction between many 



neutrinos, we study, following plR , the Schrodinger equation of a system of many 
neutrino plane waves in a box. We consider the free streaming regime and neglect 
vacuum oscillation terms, which can be included by a straightforward generalization 
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of our results. We show that the problem is equivalent to that of a system of spins, 
for which, in the case of constant spin-spin coupling, the equilibration time can be 
determined completely analytically. For neutrinos initially in flavor eigenstates, in the 
limit of the infinite number of neutrinos, this solution shows no coherent conversion, in 
agreement with our earlier results. 

The paper is organized as follows. In Sect. ^|we lay the foundations of our analysis, 
by discussing the structure of the neutrino-neutrino interaction in the flavor space and 
in the real space. We also perform an elementary analysis of coherence in a system 
of several interacting neutrinos and use it to motivate the following study. In Sect. |3| 
the many body approaches are discussed and the relevant results of our previous paper 
are reviewed. In Sect. ^ we formulate the problem of interacting neutrino plane waves 
in a box, present the full analytical solution and give its derivation. A discussion and 
conclusions follow in Sect. || and || 



2. General considerations 



2.1 The flavor structure of the interaction 



In the range of neutrino energies that are relevant to the applications (~ 10 6 — 10 7 eV), 
the interaction between neutrinos is described by the low-energy neutral current (NC) 
Hamiltonian 



H 



NC 




w2 v bi^b\- (2.1) 



Higher order effects, such as those coming from the expansion of the propagator or 
particle production (bremsstrahlung) will be neglected in our analysis. 

A property of this Hamiltonian that will play a very important role later is its 
invariance under the rotation of the SU{2) flavor group. Because of this property, the 
interaction between pairs of neutrinos - viewed in flavor space - must be equivalent to 
the interaction between pairs of spins. 

The one-particle equations of |2j, 0, |], ^TJ indeed preserve the SU{2) structure of 



the problem and have a form of the interaction between spins 1 . It would be logically 
inconsistent, however, to rely on this fact in the analysis, in which these equations are 
being tested. Therefore, below we present an explicit proof of the equivalence, valid 
regardless of whether the interaction is coherent or not. 



1 In [^2|, this property of the equations has been used to give an elegant physical analogy designed 
to explain the "synchronized" oscillations of neutrinos in a neutrino background. 
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Let us consider the interaction of two neutrinos, and, in the interests of clarity, 
for a moment suppress the Dirac indices on the fermions and the gamma matrices. 
Since the NC interactions conserve flavor, in the Hamiltonian ( |2.1| ) the flavor space 
wavefunction of a given outgoing neutrino is equal to that of one of the incoming 
neutrinos. Let us denote the flavor wavefunctions by and $j. There are then two 
possible combinations: 

#*<%%$*4^, (2.2) 
VtSuQ&tS^. (2.3) 



Since Yli=i l^«| 2 = Yli=i 1^*1 2 = ^ ne t erm m Eq. ( |2.2| ) equals 1. The term in 
Eq. ( |2.3| ) could be transformed using the well-known property of the (4-dimensional) 
a-matrices, cr^o^f = — 25^5^ (in the notation convention of 123(1), or 



2Su8 jk = SijS k i + &ij ■ a ki - (2.4) 

One gets 

^*S a ^*5 jk ^ k = \{l + ■ n^i)- (2.5) 

In this form, the equivalence between a system of neutrinos and a system of spins is 
manifest. The complete (including both the contributions from ( |2.5|) and (|2.2|) ) flavor 
space Hamiltonian for the interaction of two neutrinos, 1 and 2, is proportional to 

(3/2 + 01 -ay 2), (2.6) 

which coincides with the square of the operator of the total angular momentum, 1? = 
(c?i/2 + (X2/2) 2 , as expected. 

The strength of the interaction depends on the scattering angle in real space. Let 
us denote the states of the incoming neutrinos as \a and ipp and those of the outgoing 
states as ^ and x'a- The spatial dependence of the scattering amplitude is given by 

f^ m ^xXXa = 1 - (?<?V) " (X'<?X) (2-7) 

If the neutrino wavepackets are sufficiently broad so that several wavepackets over- 
lap, the preceding argument can be trivially generalized to show that the interactions 
between the "overlapping" neutrinos have the structure of the interactions between the 
corresponding spins. In the case of neutrino plane waves in a box, each spin interacts 
with all the others and the total Hamiltonian is the sum of the interactions of all pairs 
(see Sect. |~J). 
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2.2 Coherence of neutrino-neutrino scattering: basic features 

We next discuss general properties of coherent neutrino-neutrino scattering. 

By definition, scattering is coherent when the waves scattered by different parti- 
cles in a target interfere with each other. When the scatterer is distinguishable from 
the incident particle, the coherence condition is satisfied when the incident particle is 
scattered forward. For other directions, scattering is incoherent (unless the particle 
spacings in the target satisfy particular conditions, e.g. periodicity). In contrast, for 
neutrino-neutrino scattering, scattering can be coherent when the momentum of the 
scattered neutrino coincides with the momentum of one of the incident neutrinos. The 
corresponding values of the scattering angle (3 are and G, where denotes the angle 
between the two initial-state neutrinos. 

For the purpose of studying coherent effects, one may replace the full Hamiltonian 
(|2.1|) with a "toy" Hamiltonian which restricts scattering angles to values and G only, 
for which coherent interference may occur. This simplified interaction can be viewed as 
a flavor exchange between neutrinos which do not change their momenta. We emphasize 
that, while all coherent effects in the system are fully captured in this way, most of the 
incoherent effects are left out. Hence, one must be careful while interpreting incoherent 
effects found in this framework, as will be seen later. 

Following 0, we write the reduced Hamiltonian for a pair of interacting neutrinos 
with momenta k and p in the form 



where V is the normalization volume. The value of the angular factor, (1 — cosG), 
follows from (|2.7|) (taking into account the anticommutativity of the spinors), while the 
flavor structure of the interaction is determined by Eq. ( |2.6| ). 

The condition (3 = G by itself, however, does not guarantee coherence. We can see 
that by considering the scattering of an electron neutrino on two "background" muon 
neutrinos. Suppose that the initial momenta of the muon neutrinos are orthogonal to 
that of the electron neutrino and, further, that the neutrino wave packets are sufficiently 
narrow so that the two scattering events are independent, with each given by Eq. ( [2.8|) . 
Then, as a result of the interaction, we get 



/ \Ve{k)Ve{p)) \ 

.£ We(kH(p)) _ Vzg f 

dt \v^(k)v e (p)) V 

\K(fcK(p)) / 



cos Q) 



(2 (A 
110 
110 

v 2 j 




(2.8) 




\F) = (1 + 2^)1^)1^^) +ia(|i/ M )|z/ e i/ At ) + 1^)1^^)), 



(2.9) 



where a = -y/2G F {l - cosQ)t/V. 
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Both the interaction time t and the volume V depend on the size of the wavepackets. 
Taking the wavepackets for simplicity to be spheres of size I, we can write 



t~l,V~l 3 , (2.10) 
a G F (1 -cos6)// 2 . (2.11) 

We notice that \a\ <C 1, as long as the size of the wavepackets is much greater than (100 
GeV) -1 , i.e., the neutrino energy is well below the weak scale. This justifies neglecting 
terms of order a 2 and higher in Eq. ( |2.9|) . 

Let us compute the probability P that, as a result of the interaction, the incident 
neutrino was converted into v^. Introduce the tl v^ number" operator p, = that 
acts only on the state of the incident neutrino. Using the orthogonality of the states 
\ueUfj) and |z^v e ), the probability in question is 

P = (F\p,\F) = 2a 2 . (2.12) 

This shows that the process is incoherent: the probability we found coincides with the 
sum of the conversion probabilities for each scattering (a coherent process would yield 
P = (2a) 2 ). As can be easily confirmed, as the number of "background" neutrinos N 
is increased, the probability of flavor conversion goes as Na 2 . 

It can be further shown that if the background neutrinos are in a flavor superposi- 
tion state u x = cosaz/ e + sinaz^, the probability of conversion is not, in general, equal 
to the sum of the probabilities for each scattering. Indeed, a straightforward calcula- 
tion yields P = a 2 sin 2 a{2 + 2 cos 2 a). Comparing this with the conversion probability 
for a single scattering event, a 2 sin 2 a, we see that the conversion process is partially 
coherent, and the degree of coherence increases with decreasing a. A generalization to 



a larger number of scatterers is EC 



P = a 2 sin 2 a((iV 2 -A0cos 2 a + A0, (2.13) 

which clearly shows the presence of partial coherence (the terms proportional to N 2 ). 

The above example shows that whether a neutrino undergoes coherent conversion 
or not depends on the flavor state of the many background neutrinos it interacts with. 
Moreover, as can be seen from Eq. ( |2.9| ), each neutrino- neutrino interaction creates an 
entangled state. Hence, the flavor evolution of a system of several interacting neutrinos 
demands a many-particle description. The question is: does the wavefunction of the 
system somehow factorize into individual one-particle wavefunctions in the limit of 
large N , and, if yes, under what conditions? Below we describe two recently proposed 
approaches to this problem. 
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3. Single- vs. many-particle description: two approaches 



3.1 First approach: interference of many elementary scattering amplitudes 

The first method was developed in [2(|. The idea is to describe the flavor evolution in a 
neutrino gas by generalizing the procedure of adding elementary scattering amplitudes, 
as was outlined in the three-neutrino example of Sect. [2.2| . 

Let us consider the interactions between two orthogonal neutrino beams. Once the 
results for this setup are understood, they can be straightforwardly generalized to more 
complicated systems. Let us assume that the neutrino wavepacket size is much smaller 

I/O 

than the particle spacing (/ <C n v ), so that different neutrino- neutrino interactions 
can be treated independently using Eq. ( |2.8|) . As an initial configuration, we take the 
first beam to be made of N\ electron neutrinos z/ e , and the second beam to be made of 
N2 neutrinos in the flavor superposition state v x = cosaz/ e + sinaz/^. It can be shown 
[ pOf that the result of the interaction during a small time 5t can be written as 

\eee...)\xxx...) *^±?* |_p^ = \eee...)\xxx...) +ia\F 1 ), (3.1) 
\Fi) = iViiV 2 (l + |(e|x)| 2 )|eee...)|xxx...) 

+ N 2 (n\x)(x\e)(\fiee...) + \efie...) + ...)\xxx...) 
+ Ni(e\x) (y\e)\eee...)(\yxx...) + \xyx...) + ...) 
+ (ii\x){y\e)(\iiee...) + \efie...) + ...){\yxx...) + \xyx...) + ...), (3.2) 

a direct generalization of Eq. Q2.9| ). Here the state \y) is defined to be orthogonal to 
\x): (x\y) = 0. 

The evolution in Eqs. ( |3.1| , |3.2| ) cannot be described using only one-particle equa- 
tions, since the final state in Eq. (|3.2j ) is an entangled many-particle state. Nevertheless, 
as observed in p0| , the coherent part of the evolution can be. Indeed, the last term in 
Eq. (|3.2|) represents an incoherent effect, since it contains the sums of mutually orthog- 
onal terms, and could be dropped. The sum of the remaining three terms is equal, to 
the first order in a, to a product of rotated single-particle states. 

Since for small 5t a state of each neutrino undergoes a small rotation and since noth- 
ing in this argument depends on the particular choice of the initial states, Eqs. ( |3.1| . |3.2| ) 
specify how the states will rotate at any point in time. This observation means that, 
for each neutrino, one can write a differential equation describing its evolution for any 
— not necessarily small — t. This equation takes on the form 

.d\^ {i) ) 



(jf H aa \^), 
H aa = V2G F n u \l - cos 6 fo) ) 



|0W)(0^| + i-i|(0W|^))| 2 



(3.3) 
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where, for generality, we restored the angular factor and summed over possible angles 
QM. The flavor -space wavefunction of the given neutrino is denoted by ip^\ those 
of the other ("background") neutrinos by cp^\ The quantity denotes the number 
density of neutrinos whose momenta make an ang le 9 fo) with the given neutrino. This 
equation is similar to the one-particle evolution equations given in [|I|, ||. In fact, the 
only difference is in the last two terms in the square brackets, which have no effect on 
flavor evolution. 

Let us discuss a particular case when the neutrinos at t = are in the flavor 
eigenstates. We note that, according to Eq. ( |3.3| ), a coherent conversion effect in this 
case is absent. This result in turn could be traced to Eq. (pT2|), where the second and 
the third terms vanish for v x = or v x = v e . The only remaining source of flavor 
transformation is the last term in Eq. Q3.2|) , which is nonzero for u x = v^. Since this 
term represents incoherent scattering, for a neutrino traveling through a gas of neutrinos 
of opposite flavor the conversion probability for small t has the form P(t) = N(t)a 2 . 

At first, it is not obvious how this incoherent effect is related to the incoherent 
flavor exchange that occurs in a real neutrino ensemble. Indeed, as mentioned in the 
introduction to Eq. (|2.8|) , by replacing the true interaction Eq. ( p.l|) with the forward 
scattering Hamiltonian, we have truncated most of the incoherent effects. Nevertheless, 
a connection between the two can be made. The argument, which follows below, is very 
instructive and will prove helpful in the subsequent Sections. 

When discussing the "forward scattering probability", we in fact mean the prob- 
ability for the neutrino to be scattered in a certain (small) solid angle SQ around the 
direction of the incoming neutrino (the probability of neutrino scattering at any fixed 
angle is, of course, zero). A neutrino wave packet with the cross section area A and 
wavelength A = 2ir/E v is diverging due to diffraction into a solid angle SQ^is = \ 2 / A. 
The neutrino wave which is scattered into the same solid angle SQ will interfere with 
the incident wave. Thus, 5Q = SQdm- 

With this in mind, let us consider the result of forward scattering of a given neutrino 
on N other neutrinos in the medium. According to Eq. (|2.8| ), for neutrino wave packets 
of longitudinal size I and cross section A, a ~ — G pi (1 — cos Q)/ Al = — Gp(l — cos0)/A 
The resulting flavor conversion probability P depends on whether the forward scattering 
amplitudes add up coherently or incoherently. In the first case, P co h ~ N 2 a 2 , in the 
second, P; nc ~ Na 2 . The number of interactions during time t is iV = Atn u [n v is the 
neutrino number density) and we find that in the case of coherent forward scattering 
the dependence on the size of the wave packet drops out, 

Pcoh ~ G|(1 ~ 2 C ° Se)2 ^V = G 2 F (l ~ cos Q) 2 t 2 nl. (3.4) 
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This is indeed seen from Eq. ( |3.3|) , which does not contain the size of the neutrino 
wavepacket. 

On the other hand, for the case of incoherent forward scattering, the cross section 
of the wave packet does enter the final result, 

G 2 (l-cos6) 2 A G 2 (l -cos6) 2 
^inc ~ FK A2 -Atn v = M Un u . (3.5) 

The dependence on the size of the wavepacket look puzzling. Nevertheless, it is easily 
understood once we recall that, by construction, the result in ( |3.5| ) represents the 
fraction of incoherent scattering events for which one of the two neutrinos is scattered in 
the forward cone SQ ~ X 2 / A oc 1/E 2 A. This means that, up to a numerical coefficient, 
the probability of incoherent scattering in any direction P™ T is given by 

i^ c OT ~ Pinc/SQ ~ G 2 F E 2 U {1 - cosQ) 2 tn u . (3.6) 

This has the form a w tn u , where a w is a typical weak interaction cross section, a w oc 
G 2 F E 2 , exactly as one would expect. 

Thus, in summary, for the purpose of studying coherent effects, the full Hamiltonian 
(|2.1|) may be replaced by the forward scattering flavor exchange Hamiltonian, as done 
in Eq. Q2.8|) , with the understanding that in the final result the terms that go like N 2 
represent the complete coherent effect while the terms that go like N correspond to only 
the fraction of the incoherent scattering that occurs in the forward diffraction cone. 

3.2 Second approach: solving a many-particle equation for neutrino plane 
waves in a box 

An alternative method of studying coherent processes in a neutrino gas was proposed in 
IfHj . The idea is to consider a system of interacting neutrino plane waves in a box and 
study its flavor evolution, again with the restriction to forward scattering only. The 
neutrinos are taken to be initially in flavor states and the goal is to test the existence 
of coherent conversion effects by looking at the rate of equilibration of the system. 

The "forward scattering" Hamiltonian can be found by generalizing Eq. ( |2.8| ). A 
basis of states is formed by the initial configuration of neutrinos and all its possible per- 
mutations. Taking, e.g., N electron and N muon neutrinos, there are (2N)\/(N\) 2 dis- 
tinct basis states. The entries of the Hamiltonian are found as follows, (i) The diagonal 
entries receive N(3N-2) contributions (this includes (2iV)!/[2!(2iV-2)!] = N(2N-1) 
flavor-blind scattering processes with scattering angle (3 = and 2N\/[2\(N — 2)!] = 
N(N — 1) processes with scattering angle (3 = O). (ii) Each off-diagonal entry receives a 
single contribution from the exchange process that connects the two corresponding basis 
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states; if the two basis states are not connected by a single permutation, the correspond- 
ing off-diagonal entry vanishes, (iii) Each contribution equals \/2Gf(^ — cosOj^/V, 
where 0^ is the angle between the momenta of the two interacting neutrinos. 

For example, for the case of N = 2 there are 6 basis states: li^ei^v^}, l^z^e^); 
K^e^e), WeWeVfj), We^^e) , W^^e^e), and the Hamiltonian is given by 





^12,34 
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fu 
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^12,34 
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(3.7) 



Here = 1 — cosOj.,- is the angular factor for the i,j pair and dy t jy = + + J^, 
where is the sum of all distinct 



Let us discuss the meaning of the equilibration time i e q- If conversion is coherent, 
t eq is defined as the time over which a neutrino of a given momentum is expected 
to change its flavor. The momentum itself stays unchanged over this time scale. If 
conversion is incoherent, t eq is defined as the time over which the neutrino momenta 
will be randomized (each neutrino can no longer be identified by its momentum). In the 
model we are considering, however, the definition of equilibration for coherent process 
also applies to the incoherent process, since the interaction is chosen to preserve the 
neutrino momentum. 

Let us determine what dependence of t eq would be a sign of coherent or incoherent 
conversion. To do this, it is useful to repeat the exercise of Sect. [3.1| of estimating the 
conversion probability as a function of time. For neutrino plane waves in the box the 
following two modifications to the argument need to be made: first, the normalization 
volume V becomes the size of the box and, second, the interaction time becomes the 
actual time t elapsed from the beginning of the evolution. The scattering amplitude 
becomes a ~ — Gpt{l — cos Q)/V and, since each neutrino simultaneously interact with 
all the others, the number of interactions is iV = Vn v . Using these ingredients, we find 

P„h ~ GMl " y C ° S9>2f V - n ii = G|(l - cos e)V-nl (3.8) 

Ptec „ C,(1-«.9)V = ^(l- co.8)' 
v A V 

Here, as in the case of the interacting neutrino wavepackets, we find that in P co h the 
dependence on the size of the box cancels out (and the result is the same as Eq. (|3 . 4j) ) . 
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while for incoherent scattering the size of the box enters the final result. Since, as was 
already stressed, the incoherent scattering result is, in some sense, artificial, there is no 
cause for alarm. 

Inverting Eqs. ( |3.8| , |3.9| ) and Eqs. ( |3.4| , |3.5| ), we find that the equilibration time for 
coherent scattering in both cases has the dependence 

t c e ° q h = {G F n v y\ (3.10) 

while for incoherent scattering one finds 

for the case of neutrino wavepackets, and 

q c = (GVv/^AT 1 (3.12) 

for plane waves in a box. In what follows, we solve the evolution of the system of 
neutrinos in the box completely analytically. This gives the equilibration time t eq , 
which we compare to Eqs. (|3.10|) and (|3.12|) to determine the character of the evolution. 



4. Plane waves in a box: analytical treatment 
4.1 Setup 

As done in ||21|1 , we consider a system of 2N neutrino plane waves, of which N are 



initially in the v e state and the remaining N are in the state. We write this config- 
uration of the system as: 

|5(0)) = | u e v e ...v e i>„u ti ...P IJ ) . (4.1) 



Here each neutrino wave is identified by its position in the list, which corresponds to a 
given momentum. 



The Hamiltonian for this system has the form described in Sect. [3.2| . Since it is 
equivalent to the Hamiltonian of a system of spins (Sect. |2.1|) , in the following we adopt 
the terminology of angular momenta. As a convention, we identify the "up" (pointing 
along the z axis) state of a spin with v e and the "down" state with v^. Thus, the initial 
state ([O]) can be expressed as 

15(0)) = | + i, + 1,-1, -!,...-!). (4.2) 

S v ' V v ' 

N N 
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To make an analytical treatment possible, we set 



h = 1 (4-3) 

for all neutrino pairs. The resulting Hamiltonian describes an ensemble of spins in- 
teracting with each other with equal strength 2 . This simplifying assumption does not 
change the coherent (or incoherent) character of the conversion effects, as also pointed 
out in ||21|| . With this assumption, the nonzero off-diagonal entries of the Hamiltonian 



become equal to \/2Gf/V while the diagonal entries become V2Gf/V x N(3N — 2). 

We recall (Sect. |3.1| ) that for the initial configuration (|4.1|) the single-particle 
formalism predicts no coherent effects, as can be seen from Eq. ( |3.3| ). We investigate 
if this conclusion is changed once the simultaneous interactions between all neutrinos 
are considered. 

4.2 Results 

Let us first display our results. Their derivation and a more detailed discussion is given 
in Sec. [L3|. 

The probability Pi that the first spin - which is in the state | + 1/2) at t = 
- is found in the same up state at time t (t > 0) can be calculated by a completely 
analytical procedure. It is given by the Fourier series: 



1 1 N ~ X 

AW = - + -^ V(N, J) cos[(2J + 2)gt] , (4.4) 



j=o 



where g = \f2GpjV . The coefficients r](N, J) are the squares of the Clebsch-Gordan 
coefficients (jimij 2 m 2 1 JM) = C(ji,mi, j 2 ,m 2 , J, M) for the addition of two angular 
momenta, j\ = m,\ = (N — l)/2 and j 2 = —Tn 2 = N/2: 



r)(N, J) 



' fN — I N-l N N j 1 l x 12 

V 2 ' 2 '2~'~2~' + 2'~2 
' (l + 2J)iV[(iV-l)!] 2 (4 _ 5) 



(J + N+ 1)!(JV- J- 1)!' 



We follow the standard notation, in which ji,j 2 denote the values of the angular mo- 
menta being added, mi,m 2 denote their projections along the z direction, and J and 
M are the corresponding values for the total spin. 

Figure |l| shows the survival probability Pi(t), obtained by summing numerically 
the series fl4.4j). The horizontal axis represents time in units of the equilibration time 



2 This makes it different from the familiar Ising model. 
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Figure 1: The time-evolution of the survival probability Pi, for different values of ./V (num- 
bers on the curves). The unit of the horizontal axis is chosen to be tjt c °q = 2tgN, as used in 
Iffij] . If the conversion were coherent, the curves should reach equilibration at the same point 
on this scale. 



predicted for coherent conversion, t c ° h = (2Ng)~ 1 (Eq. ( [3T0D for 2N neutrinos). The 



curves refer to iV = 40, 160, 640, 2560. Having the same choice of units, our plot can 
be directly compared to fig. 2 of ref. plf , where, however, smaller values of N were 
used 3 . From Fig. [I], it is clear that the equilibration time t eq is larger than t c °^ and the 
discrepancy grows with N. 

In the limit N ^> 1 the series ( (4.4|) is well approximated by the following function: 



(4.6) 



P x {t) = 1 - ^VNgtexp(-Ng 2 t 2 )erti(VNgt), 



(4.7) 



where erfi(z) is the imaginary error function 

2 f z 

erfi(,2) = — i eii{iz) = —= x / exp(t )dt . 
Jo 

3 In ref. pl[ | a different form of the spin-spin (neutrino-neutrino) coupling was adopted, namely, 
the angular factors were randomly generated. This produces a numerical difference with respect to 
our results, but does not affect the conclusions about the equilibration time. 
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We have checked that for the values of N shown in Fig. [I] the approximation is very 
good and the curves obtained using (|4.6| ) are indistinguishable from those obtained 
from Eq. (p~4"l). 

From Eq. ([4.6|) we see that the characteristic equilibration time equals: 



(gVN)- 1 . (4i 



This scaling of the equilibration time can be clearly seen in Fig. [1]. It is precisely (up to 
numerical factors) the form ( 3.1 2|) expected if the equilibration process is incoherent. 



4.3 Derivation of results 

4.3.1 The eigenvalues of the Hamiltonian 

From Eqs. ( p.6|) and ( |4.3|) it follows that for each pair i and j the interaction is given 
by Hij = gL^j = g{pi + ffj) 2 and, therefore, the Hamiltonian for the whole system is 

2JV-1 2N /27V-1 27V „ \ 

H N = g j2 E (** + = 9 E E + 2 N{2N ~ 1} ' (4 - 9) 

j=l j'=j+l \ i=l j'=i+l / 

where the summation is done over all N(2N — 1) pairs. 

We notice that this Hamiltonian can be related to the operator of the square of the 
total angular momentum of the system, 

(27V \ 2 2JV-1 27V 
E^O = E E 2^ + -iV. (4.10) 
i=l / i=l j=i+l 

By comparing Eqs. ( |4.9| ) and ( |4.10|) we find 

H N = g[L 2 + 3N(N-l)\, (4.11) 

which has the eigenvalues 

E NfJ = g[J(J+i)+3N(N-l)], with J = 0, 1, iV . (4.12) 

Since the flavor evolution will depend only on the difference of the eigenvalues (as will 
be shown later), in the following we will omit the second term in Eq. ( |4.12j ) and write 

Ej = gJ(J+l) . (4.13) 

A given value of Ej in general corresponds to more than one state. This degeneracy 
can be seen by simply counting the states. For our specific setup, with M = 0, there are 
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(2N)\/(N\) 2 basis states, while there are only N + 1 different values of J. In general, 
therefore, an orthonormal basis of eigenstates of the Hamiltonian has the form: 



(4.14) 



where the index kj labels eigenstates with the same value of J. 

As an illustration of the present discussion, we give the Hamiltonian for the case 
N = 2 of Sect. 

/8 1 1 1 1 0\ 
18 110 1 
118 11 
110 8 11 
10 118 1 
\Q 1 1 1 1 8/ 



H 



N=2 



(4.15) 



The eigenvalues, which are degenerate as expected, are #(6,6,8,8,8,12). They corre- 
spond to the following values of J: (0, 0, 1, 1, 1, 2), as can be seen from Eq. ( |4.12|) . 

4.3.2 The Fourier series 

Let us start with proving Eq. ( f4.4|) . 

In the interest of clarity, we first display some definitions. Let us denote as | + 1/2) i 
and | — l/2)i the "up" and "down" states of the first spin. Given a state \ip) of the 
2N spins of the system, it will be convenient to use the projections (+l/2| 1 |?/>) and 
(— l/2|i| , 0), which are states of 2N — 1 spins. In particular, we introduce the state \S~) 
obtained by removing the first spin from the initial configuration ^(O)) (Eq. f|L2] 



(+1/21x15(0)) = \S-) 



1 

*2' 



1 

'2' 



1 



1 

2' 



I -I) 
2' 2 ; ' 



(4.16) 



— v — 

A 1 



-v- 

A' 



For any state, we can define another state obtained from the first one upon a 
reflection along the z-axis. In this transformation, each constituent spin | + 1/2) is 
transformed into spin | — 1/2) and vice versa. For instance, the application of this 
transformation to \S~) gives its "mirror" state: 



\S + ) 



1 

2' 



1 



1 l ^ 

2- + 2 ) 



(4.17) 



N-l 



N 



In what follows we will exploit the properties of various states under this reflection. 
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We adopt the notation \j,m, s) n to indicate the result of the following procedure: 
consider the state \s) and take its projection on the subspace formed by the states 
of n spins that have the total angular momentum j and projection m along the z 
axis; normalize the result to 1. As an example, |J, 0, S(0))2N indicates the normalized 
projection of 5*1(0)) on the subspace of the states of 2N spins having total angular 
momentum J and ^-projection equal to zero. 

It is also convenient to define the state: 

|sj> = <+l/2| 1 |J,0,5(0)) 2JV , (4.18) 
(which has 2N — 1 spins) and denote as \sj) its mirror state obtained upon ^-reflection. 

The proof is constituted by four parts. 

1. Expansion of the initial state 



The initial state can be expanded on the basis of the eigenstates of the Hamiltonian: 

\S(0)) = Y,\Wj)(J,kj\S(0)h ( 4 -!9) 

J,kj 

which then can be evolved in time as 

\S(t)) = Y,e- lEj V,kj)(J,kj\S(0)). (4.20) 

J,kj 

Since the energy depends only on J, the summation over kj can be done, yielding 

N 

|5(0)> = J2 A Mo,s(o)) 2N , 

J=0 

N 

\S(t)) = J2^je- lEj V,0,S(0))2N 

The factors Aj equal the Clebsch-Gordan coefficients: 

„ (N N N N \ 

This can be understood considering that the initial state, 15(0)), equals the sum of two 
angular momenta with j 1 = N/2, j 2 = N/2 and projections m\ = N/2, m 2 = —N/2. 
The first momentum is given by the sum of all the "up" spins, while the second is the 



(4.21) 



(4.22) 
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sum of all the "down" spins. 

2. Decomposition into states of 1 and 2iV — 1 spins 



It is useful to decompose the state \S ) as: 



N-l 

\S~) = Kj+iMJ+ "2' S ~)™-i • ( 4 ' 23 ) 



j=o 

Similarly to Eq. (|4.22|) , here we have: 



Km , = c (!Lz±, !Lz±, f , -f , j + i -i ) = ^7) , (4.24) 



where rj(N,J) was given in Eq. ( |4.5| ). The argument to explain ( 4.24 ) is analogous 
to that used to justify Eq. ( |4.22[ ): since the state \S~) has N — 1 spins "up" and N 
spins "down", it is given by the sum of two angular momenta with ji = (N — l)/2, 
j 2 = N/2 and projections mi = (N — 1)/2, r«2 = —N/2. These momenta are the result 
of summing all the spins "up" and all the spins "down" separately. 
Each of the states \J, 0, S(0))2N can be decomposed as: 

\J,0,S(0)) 2 n = aj\ + -)i|J + -,--,Sj) 2N -i + h\ + — >x|«^ — -,--,s7) 2 iv-i 

+ cj\ - -)x| J + i +i 4) 27V -i + dj| - J - ^, 4>2iv-i (4.25) 

where the states |J + |, +33 Sj) 2 n-i and |J — |, +|, Sj) 2 at_i in the second line are 
related to the corresponding ones in the first line by the ^-reflection operation. Notice 
that the state | J, 0, S(0)) 2 n has parity (— 1) N ~ J under ^-reflection. Requiring that the 
right hand side of Eq. (|4.25 ) has the same parity yields 



CJ = (-1) N - J aj dj = (-1) N - J bj . (4.26) 

We also note that 

a N = c N = 6 = d = , (4.27) 

\aj\ 2 + \bj\ 2 =\cj\ 2 + \dj\ 2 = 1/2, (4.28) 

where the last relation represents the normalization condition. 
3. Properties of the decomposition ( 4.25 ). 
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We can calculate the state (— l/2|i|S(£)) and demand that it vanishes at t — 0. 
This is equivalent to requiring a null probability to find the first spin initially in the 
down state. Combining this condition with Eqs. (|4.21| ), (|4.25| ) and ( |4.26| ) we find 

N-l 

j=o 



N-l 



J+1<JJ+1\ 



1) 



N-J-l 



J=0 



1 1 

J + 2' + 2' S J+^ 2N ~ X = ' 



from which it follows that 



Ajaj\J + -, +~, s)) 2 tv-i = A J+1 b J+1 \J + ^, +~, s^ +1 ) 2 n-i , 



and, by ^-reflection: 



1 



Ajaj\J + -,--,Sj) 



2N-1 



2' 2 

4. Calculation of the survival probability Pi. 



Aj + ib J+1 \J + -, --, s j+1 ) 2 at-i 



(4.29) 



(4.30) 



(4.31) 



Let us calculate the probability -Pi(t) that at the time t the first spin is found in 
its initial state, | + 1/2): 

P 1 (t) = (S(t)\ + l/2) l (+l/2\ 1 \S(t)) . (4.32) 
The combination of Eqs. (|PID , ( ]4T25| ) and ( P~2TD , gives the expression of (+l/2|i|5(t)): 

N-l 



(+1/Wt)) = e- tEjt Ajaj\J+ I -I, S j) 2N ^ 

N-l 1 l 

+ e- lE ^ t A J+1 bj +1 \J+ -, --, Sj +1 ) 2N ^ 



(4.33) 



j=o 



which, taking into account the result ( f4.31|) , simplifies to: 

N-l 



1 1 



j=o 



(4.34) 



The final step to conclude the proof is the observation that for t = 0, Eq. (|4.34|) 
should coincide with the expansion of the state \S~), given in Eq. ( [4.23|) . From 
comparison, using the equality ( |4.24|) , it follows that: 



1 1 1 1 

2Ajaj\ J + -, --, Sj) 2 n-i = K J+ i/ 2 \J + -, --, S ) 2 n-i 

= ^N^jj\J + ±-l,S-) 2N -i 



(4.35) 
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Combining Eqs. flO%) , flOl , (|435|) and (PD , we find 



= ^ ^(^ J ) 2 | e -^ + e -^ +1 t| 2) (4 36) 



4 
j=o 



from which Eq. ( |4.4| ) straightforwardly follows. 

4.3.3 The physical (large N) limit: resummed function 



To derive Eq. ( f4.6|) we first elaborate r)(N,J), Eq. ( |4.5| ), in the limit N ^> 1. We 
approximate the factorials using the Stirling formula: z\ ~ v / 27re (2: ~ 1 / 2 ' ) lnz ~ z . As a 



result, an exponential factor appears in the expression of r](N,J). We expand this 
exponential in powers of J/N and keep only the leading terms. This is justified because 
the various factors in Eq. ( ^4.5| ) are significantly different from zero only for J ^ N. 
From the procedure described, we obtain: 

rj(N, J) ~ 2iV~ 1 (l + J) exp [—(J + 1) 2 /N] . (4.37) 



This function has a peak at J = a/ N/2 — 1 and has a characteristic width 2a ~ y/~N. 
This fact provides an a posteriori justification for taking J iV for large N. 

As a second step, we notice that, for N ^> 1 the Fourier sum in Eq. (O) is well 
approximated by the integral: 



N-l 



^ r/(iV, J) cos[(2J + 2)gt] ~ / dJ — Je~ j2/N cos(2Jgt) 

7— a J 

= 1 - v / ^iV^te- 7V92 * 2 erfi(v / Ngt), (4.38) 



where the expression (|4.37| ) for rj has been used. The function erfi(z) is given in Eq. 



(O). From Eq. ([4.38|) and (O) the final result (KB) follows immediately. 



5. Discussion 



1. Let us compare the main features of our solution, Eqs. ( f4.4| ) and (|4.6j), with the 



corresponding results for a gas of small neutrino wavepackets. In the latter case the 
survival probability exhibits a different dependence on time, characterized by a simple 
exponential form with the argument linear in t. To see this, consider a thought ex- 
periment in which a beam of 2N electron neutrinos enters the medium containing an 
equal mixture of electron and muon neutrinos. It is straightforward to see that, for 
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incoherent forward scattering, the flavor composition of the beam will change according 
to 



dN e /dt 
dN^/dt 



A(iV e -i\g, 
A(iV M -iV e ), 



(5.1) 
(5.2) 



where A denotes the frequency of the flavor changing collisions for a given neutrino. 
The solution is 



a simple exponential. The origin of the difference between Eq. (|4.6| ) and (|5.3|) lies in 
the different physics of the interaction. In the case of the wavepackets, the time elapsed 
from t = tells us how many background neutrinos a given neutrino interacted with; 
the interaction time between neutrino pairs is very short. On the other hand, in the 
case of the interacting plane waves, any two waves in the box continuously interact with 
each other starting from t = 0. The scattering amplitude (and the effective interaction 
cross section) grows with time. If the conversion process were coherent, the survival 
probability would not depend on the particular model adopted (see Sect. [OP . 

2. It is interesting to note that the Fourier series in Eq. fl4.4j) describes a periodic 
function with period T = n/g. This periodicity is lost when a transition to the integral 
is made. Indeed, in Eq. ( |4.38|) the lowest frequency over which the integration is per- 
formed is zero (J = 0), corresponding to infinite period. It is clear that the period T 
is much larger than the equilibration time t eq (see Eq. ( ^4.8| )): T = n\/Nt eq . Moreover, 
one can see that typically T exceeds the age of the universe: for a volume V = 1 cm 3 
we find T ~ 10 22 s. Therefore the integral form ( |4.(j| ) is an accurate approximation of 
the exact result over time scales (at least) of the order of t eq , and applies to all the 
realistic physical situations. 

3. Let us comment on the range of validity of our results and their possible exten- 
sions. 

(a) Due to the SU(2) invariance of the problem, the results are the same if we 
replace the flavor eigenstates u e , with any pair of states in the flavor space, v x ,v y ^ 
related to the flavor basis by an SU(2) transformation. Explicitly, if we take neutrinos 
initially in the states v x and z/ y , the probability of conversion u x — > v y does not exhibit 
coherent effects in the large N limit. This tells us that no coherent conversion is realized 
if coherence is initially absent in the system, i.e. if, for a given neutrino Vi being in the 
flavor state v x , the other neutrinos Uj in the ensemble are not coherent superpositions 
of v x and v y : {vi\fj) = or |(z^|^j)| = 1. 



N e (t) = N(l + exp(-Xt)), 



(5.3) 



21 



(b) Several generalizations of our calculation are possible. For instance, one could 
investigate how the results depend on the initial configuration of the system. In pO 



we 



have studied the case in which coherence is initially present in the system (|(z/j|z/j)| < 1) 
and/or the initial state is entangled. The findings were in agreement with the one- 
particle description, Eq. (|3.3| ). 

(c) The assumption of constant spin-spin coupling, = 1, could be relaxed to 
describe the case of, e.g., randomly distributed couplings or interactions which de- 
pend on the spin-spin distance or on the specific geometry of the system. These cases 
may be difficult or impossible to describe analytically; however we expect them to be 
characterized by the same incoherent character of the conversion. 

(d) The generalization to include vacuum oscillations and interaction with ordinary 
matter can done by adding the appropriate (well known) terms to the Hamiltonian. 

6. Conclusions 

We have found that a system of many (N ^> 1) (forward-)interacting neutrino plane 
waves in a box reaches flavor equilibration after a time t ~ t eq — (\^2GpVN /V)~ l if 
the neutrinos are initially in flavor states. This result agrees with what is expected from 
incoherent scattering for this system, and therefore indicates the absence of coherent 
conversion effects. 

Although our results were obtained under specific simplifying conditions (plane- 
waves in a box, constant couplings, forward scattering only, etc.), they can be used to 
draw important general conclusions on neutrino conversion in a neutrino background. 
The absence of coherent conversion has general character, since it does not depend on 
the size of the neutrino wavepackets (see Sect. |3"H] ) . It contrasts with ref. []2T| and agrees 



with the prediction of the traditional - and widely applied - single-particle description. 
This tells us that there is no reason to believe in a breakdown of this description and 
that, therefore, its applications have to be considered valid. The significance of this 
statement is clear in consideration of the important physical problems, in cosmology 
and in stellar physics, in which the neutrino-neutrino coherent scattering plays a role. 

Finally, we remark that, in principle, the analytical solution we have found for a 
system of interacting spins may be interesting on its own and may find applications 
beyond the neutrino evolution problem considered here. 
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